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Using the disconjugacy properties of the Schrodinger equation, it is possible to develop a new 
type of generalized SUSY QM partnership which allows to generate new solvable rational extensions 
for translationally shape invariant potentials having a finite bound state spectrum. 

For this we prolong the dispersion relation relating the energy to the quantum number out of the 

physical domain until a disconjugacy sector. The prolonged excited states Riccati-Schrodinger (RS) 

functions are used to build Darboux-Backlund transforms which give regular isospectral extensions of 

*vj ", the initial potential. We give the spectra of these extensions in terms of new orthogonal polynomials 

I ■ and study their shape invariance properties. 
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I. INTRODUCTION 



The last three years have seen a substantial development of research works concerning the study of rational ex- 

j— r tensions of solvable quantum potentials, in particular because of their intimate relation with the recently discovered 

^G ' exceptional orthogonal polynomials (EOP) |4427j|. 
Mh, In [244271 we have proposed a new scheme to generate the rational extensions of every primary translationally 

shape- invariant potentials (PTSIP) [l|-y, 12a]. They are obtained via Darboux-Backlund Transformations (DBT) 
based on negative eigenfunctions built from excited states of the initial hamiltonian. The regularity of corresponding 
Riccati-Schrodinger (RS) functions [28j is then directly verified by combining disconjugacy theorems and asymptotic 
analysis of the eigenfunctions. This approach which is systematic and generalizes the usual SUSY QM partnership 
can also be enlarged in a multistep version. In (2J-[27J these unphysical eigenfunctions are deduced from the excited 

,_^ ' states via discrete symmetries acting on the set of parameters of the initial potential. For the isotonic or trigonometric 
£> , Darboux-Poschl- Teller potentials, when the obtained extensions are strictly isospectral to the original potential, their 

OS ' eigenstates are (up to a gauge factor) the EOP discovered by Gomez- Ullate, Kamran and Milson |4(. 

For potentials presenting an infinite number of bound states associated to all integer values of the quantum number, 
the recourse to discrete symmetries is necessary to reach the disconjugacy sectors. However for PTSIP with a finite 
bound state spectrum {E n , n < n max } there exists another way. In this case, some disconjugacy sectors can be 

C*~) . attained by prolonging the eigenstates ipn for values of the quantum number n going beyond n max . If the "dispersion 
relation" , ie E n as a function of n, goes then to negative values, the corresponding prolonged eigenstates ipn, although 
diverging at least at one extremity of the definition interval, may possibly be used to build regular extensions via 
DBT. Here again, disconjugacy theorems and asymptotic analysis allow to control the regularity of the corresponding 
RS function w n . In the present paper we consider the four PTSIP for which the associated dispersion relation has 
the mentioned behaviour. We can share them into two groups: the Morse and hyperbolic Darboux-Poschl- Teller 
potentials, which have a parabolic (quadratic) dispersion relation onb the one hand and the Eckart and hyperbolic 
Rosen-Morse (HRM) potential for which the dispersion relation has the form of a second degree Laurent polynomial 
on the other hand. For the first group the prongation leads to only one disconjugacy sector and the extensions 
obtained are strictly isosectral to the initial potential. For the Morse potential we recover exactly the new extended 
potentials obtained very recently by Quesne [291 ] . Quesne shows in particular that if these extensions don't inherit 
of the translational shape invariance of the original potential (which one of the feature of the extensions associated 
to the EOP), they satisfy a kind of "enlarged shape invariance property". We prove that it is also the case for the 
similar extensions of the HDPT potential. For the Eckart potential, again we have only one disconjugacy sector 
with strictly isospectral extensions. This is no more the case for the HRM potential. In this last case, we have three 
distinct disconjugacy sectors, the two first corresponding to strictly isospectral extensions while the DBT built on 
eigenfunctions of the third sector are reverse SUSY QM partnership and then give only quasi isospectral extensions. 
For all the obtained strictly isospectral extensions, we give explicit expressions for the eigenstates in terms of new 
orthogonal polynomials. 



II. DISCONJUGACY AND REGULAR EXTENSIONS OF ONE DIMENSIONAL POTENTIALS 

If tp\(x; a) is an eigenstate of H{a) = —d 2 /dx 2 + V(x; a), a £ R m , x £ I C R, associated to the eigenvalue E\(a) 
(E o (a) = 0) 

ip^(x; a) + (E x (a) - V(x; a)) ijj x (x; a) = 0, (1) 

then the Riccati-Schrodinger (RS) function w\(x;a) = —ip' x (x;a)/i/)\(x;a) satisfies the corresponding Riccati- 
Schrodinger (RS) equation [28| 

- w' x (x; a) + w 2 x (x; a) = V(x; a) - E\(a). (2) 

The set of Riccati-Schrodinger equations is preserved by the Darboux-Backlund Transformations (DBT), which are 
built from any solution w u (x; a) of the initial RS equation Eq@ as [2|| (32, |33[ 



/ ,A(w) (v) E x (a)-E u (a) 

w x (x;a) -> w\'(x;a) = -w v (x;a) -\ -, (3) 

w v (x\a) — w\(x; a) 

where E\ (a) > E u (a) . w x is then a solution of the RS equation: 

-w[ v) '(x;a) + (w^(x;a)) 2 = V^(x;a) - E x (a), (4) 

with the same energy E\(a) as in Eq([2]) but with a modified potential 

V {v \x;a) = V{x;a) + 2w' u {x;a), (5) 

which is called an extension of V . 

The corresponding eigenstate of H^ u \a) = —d 2 /dx 2 + V^'(x; a) can be written 

tfW, - exp (- f^V,a)) ~ ' =%) W ,;,), ,6) 

V J J y/E\(a)- E v {a) 

where A (a) is a first order operator given by 

A (w„) = d/dx + w v {x; a). (7) 

From V, the DBT generates a new potential V^"' (quasi) isospectral to the original one and its eigenfunctions 
are directly obtained from those of V via Eq©. Nevertheless, in general, w v (x;a) and the transformed potential 
V^(x;a) are singular at the nodes of ip v (x;a). For instance, if ip n (x;a) (v = n) is a bound state of H(a), V^ is 
regular only when n = 0, that is when ipn=o is the ground state of H, and we recover the usual SUSY partnership in 
quantum mechanics. 

We can however envisage to use any other regular solution of Eq([2]) as long as it has no zero on the considered real 
interval /, even if it does not correspond to a physical state. In particular, it is possible to use some eigenfunctions 
ip v associated to negative eigenvalues E v {a) < [3J|. This is due to the disconjugacy of the Schrodinger equation for 
these eigenvalues [27| . More precisely, a second order differential equation like Eq(fT]) is said to be disconjugated on I 
if every solution of this equation has at most one zero on I [35l - l37l |. For a closed or open interval /, the disconjugacy 
of Eq([T]) is equivalent to the existence of solutions of this equation which are everywhere non zero on / [35I437J . 

We have also the following result 

Theorem [35|, [36j If there exists a continuously differentiable solution on / of the Riccati inequation 

-w'{x)+w 2 {x) + G{x) <0 (8) 



then the equation 

i//'(x) + G{x)tp{x) = 0, (9) 

is disconjugated on I. 

In our case, if E\(a) < 0, we have 

— w' (x; a) + u>o(x; a) = V(x; a) < V(x; a) — E\(a), (10) 

wq(x] a) being continuously differentiable on I. The above theorem ensures the existence of nodeless solutions 4>\(x; a) 
of Eq([T|), that is, of corresponding regular RS functions v\(x; a) — —4>' x (x; a)/4>\{x; a). To prove that a given solution 
<p\(x;a) belongs to this category, it is sufficient to determine the signs of the limit values 4>\(0 + ;a) and (f>\(+oo]a). 
If they arc identical then <j)\ is nodeless and if they are opposite (f>\ presents then a unique zero on /. In the first case 
V(x; a) + 2v' x (x; a) constitutes a perfectly regular (quasi)isospectral extension of V(x; a). 

Of particular interest is to obtain such solutions 4>\(x\ a) which are polynomials (up to a gauge factor) in order 
to build rational extensions of the initial potential V. In |25l - t27| . e have shown that it is possible to generate such 
solutions by using specific discrete symmetries I\ which are covariance transformations for the considered family of 
potentials 

{ V(x; a) ±4 V(x; a. t ) = V(x; a) + U (a) . 

Ti acts on the parameters of the potential and transforms the RS function of a physical excited eigenstate w n 
into a unphysical RS function v n ^(x;a) = Ti (w n (x;a)) — w n (x;ai) associated to the negative eigenvalue £ n .i{o) = 
r t {E n {a)) = U(a)~E n {a t )<0.' 

- v' n>i (x; a) + vl,(x; a) = V(x; a) - £ n>i (a). (12) 

To v ni corresponds an unphysical eigenfunction of H(a) 



<p n ,i(x;a) = exp ( - / dxv, hi (x;a) I (13) 

associated to the eigenvalue E n> i{a). 

Since £ n ^{a) < 0, if (j) n ^{x;a) has the same sign at both extremities of /, then v n ^(x;a) can be used to build a 
regular extended potential (see Eq© and Eq([6])) 

V^(x;a)=V(x;a)+2v' n /x;a) (14) 

(quasi)isospectral to V(x; a). The eigenstates of V^ 1 '" 1 ' are given by (see EqQ) 

{ri.i) i -, , \ E k (a)—£ ni (a) 

wl (x: a = — v„ i (x: a H r^ ) \ 

■0j™ ' (x;a) = exp (- J dxw^ 1 '''' (x;a)J — A (v n)i ) i[>k(x; a) 

for the respective energies Ek(a). 

The nature of the isospectrality depends if 1/ ' 4> n ^{x; a) satisfies or not the appropriate boundary conditions. If it 
is the case, then \/<j) n .i(x] a) is a physical eigenstate of H^ n ^{a) = —d 2 /dx 2 + V^ n,l \x;a) for the eigenvalue £ nt i(a) 
and we only have quasi-isospectrality between V (x; a) and V^ n '^ (x; a). If it is not the case, the isospectrality between 
V (n ^ (x; a) and V(x;a) is strict. 

This procedure can be viewed as a "generalized SUSY QM partner ship " where the DBT can be based on excited 
states RS functions properly regularized by the symmetry Tj. In |25l - l27| . it has been applied to exceptional PTSIP 
of the first and second categories [28j . In the particular case of the isotonic oscillator, the spectrum of the two first 



series of extensions are given (up to a gauge factor) by the exceptional Laguerre polynomials initially discovered by 
Gomez-Ullate et al jg]. 

For some potentials which have finite bound states spectrum, it is possible to use an even more direct way to 
generate negative energies eigenstates. In this case, the eigenfunctions rp n (x; a) satisfy the required Dirichlct boundary 
conditions to be acceptable physical eigenstates only for a finite number of values {0, ..., n max } of the quantum number 
n. Beyond this maximal value n max , tp n have a divergent behaviour at (at least) one extremity of the definition interval 
/. ip n has then to be rejected as eigenstate but it can still be used to build the corresponding DBT A (w n ), the energy 
E n (viewed as a function of the quantum number n) being extended to values of n greater than n ma x- If hi this 
extended domain E n reaches negative values we recover then a disconjugacy sector of the Schrodinger equation and 
ip n can be exempt of nodes. The DBT A (w n ) gives then a regular extended potential (see Eq© and Eq©) 

V (n) (x;a) = V(x;a) + 2w' n (x;a), (16) 

(quasi)isospectral to V(x;a) and its eigenstates are given by (see Eq((3J|) 

(n) / \ / \ . E k .(a) — E n (a) 

wi (x: a) = —w n (x; a) H r-M — ; ' ■. 

4> k n '(x;a) = exp (- J dxw^' (x; a)) ~ A (w n ) ipk(x;a). 

III. PRIMARY TRANSLATIONALLY SHAPE INVARIANT POTENTIALS WITH FINITE BOUD 

STATES SPECTRUM 

In [23], we have shown that all the PTSIP can be classified into two categories in which the potentials can be 
brought into a harmonic or isotonic form respectively, using a change of variable which satisfies a constant coefficient 
Riccati equation. Among them we have potentials admitting an infinite bound state spectrum [l|, |2|, |28|. This is 
naturally the case of the confining potentials as the harmonic and isotonic oscillators, the trigonometric DPT (or 
Scarf I) and the trigonometric Rosen-Morse potentials. For these ones the dispersion relation is a strictly increasing 
function on n £ [0,+oo[ which is linear for the two first, parabolic for the third one and which is a second degree 
Laurent polynomial for the fourth one. This is also the case of the effective radial Kepler- Coulomb (ERKC) potential 
which admits an infinite bound state spectrum with negative energies and a continuous spectrum of scattering states 
with positive energies. The dispersion relation for the bound states (a second degree Laurent polynomial without a 
regular term) is also strictly increasing for n£ [0, +oo[. 

The other primary TSIP have a finite bound state spectrum. This is the case in the first category for the Morse, 
Eckart and hyperbolic Rosen-Morse potentials and in the second category for the hyperbolic Darboux-Poschl- Teller 
(or Scarf II) potential [H, 0, H] . 

A. Case of dispersion relation which is a second degree polynomial 

1. Morse potential 

For the Morse potential 0, [2, [13 (x e M) 

V(y; a, b) = b 2 y 2 - 2 ( a + ]- j by + a 2 , a, b > 0, (18) 

where y = exp (— x), we have a "dispersion relation" (energy E n as a function of the quantum number n > 0) which 
is of parabolic type. Namely, we have (a n = a — n) 

E n (a) = a 2 - al = n (2a - n) . (19) 

The bound states are obtained on the increasing part of E ni that is for n < a and we have exactly [a] bound states 
( [a] being the integer part of o) given by 



iff n (x; a, b) = % fa a n , b) L 2 n a - (2by), (20) 

where 

^ (x;a,b)=y a e- b y. (21) 

The corresponding RS functions are given by 

w n (x;a,b) = w (x;a n ,b) - 2by "^ . , - , (22) 

L„ a "(2by) 

where 

wo (x; a,b) = a — by. (23) 

Beyond n — a, the ipn no longer correspond to physical eigenstates and the E n function decreases. When n exceeds 
the value 2a , E n becomes negative and the corresponding Schrodinger equation enters in a disconjugacy regime. 
Using 



x ^ 0+ i iv> . ( 24 ) 

L n \ z ) , n\ Z 



x — > + oo 



where (X) = (X) ... (X + n — 1) is the usual Pochhammcr symbol [31], we find that the asymptotic behaviour of 

tp n (x; a, b) at +oo and — oo is given by 



_ „ (2a - 2n + 1) ... (2a - n) , s 

il> n (x;a,b) ~ y n ~ ay - f— ^ i -> ±oo, (25) 

a;— >+oo 77,; x~ f+oo 



and 



Vnfoa.ft) ~ ^_i_y» e -6« _> ± , (26) 

£— > — oo 77,1 a:— >+oo 

where ± = (—1)™. 

We see that ij) n has the same sign at both ends of the definition interval, the disconjugacy of the Schrodinger 
equation that it satisfies implies that ?p n has no zero on R. We can then use the corresponding RS functions to build 
(quasi)isospectral extensions of the Morse potential respectively. They are given by 

V y ' (x; a, b) = V(y; a, b) - 2y , (27) 

ay 

where w n (y, a n , b) is given in Eq (|22[) . with y = exp(— ax). 

Since l/ij) n (x; a, b) is divergent at — oo, it is not a physical eigenstate of V^" 1 ' which is then strictly isospectral to V 
and admits for bound state eigenfunctions (fc <E {0, ..., [a — 1]}, z = 2by) 

L a (z) 
ip k n >(x;a,b) = (w n (x;a,b) - w k (x;a,b))ip k (x;a,b) = ip (x;a k ,b) "^ , (28) 

Lin n (Z) 



(29) 



where, making use of the following recurrency properties of the Laguerre polynomials [30l |31| 

L^(z) + L^(z) = L^(z) 

zL« 1 l(z) = (n + a)L«_ 1 (z)-nL? l (z), 

we can write 

L^ k (z) = (2a - k) L»- (z)Ll a _\(z) - (2a n) L 2 k a » (z)^(z), (30) 

which is a polynomial of degree n + k — 1. 



2. HDPT potential 

The dispersion relation of the hyperbolic Darboux-Poschl- Teller (HDPT) potential Q, Hill > 0) 

(g + l/2)(«-l/2) 03 + 1/2)03-1/2) , T „ ,i^i/9 rm 

K(o:;a, J 0) = — ^ ~2 hVo(a,Pj, /3>a + l>l/2, (31) 

sinn a; cosn x 

with Vb(a, /3) = (/3 — a — l) 2 , is of the same parabolic type than in the preceding case 

E n (a, 0) = An (2a - n) , (32) 

with a = (/3 — a — 1) /2. Here again the bound states are obtained on the increasing part of E n , that is, for n < a 
and they are given by 

</>„ (x; a, /?) = ^o (x; a, /3) J**'-« (z) , (33) 

where z = cosh2x and 

iPo(x;a,/3) = {smhx) a+1/2 {cosh X y f3+1/2 . (34) 

using [ad, m 

(^(^f^^^ +1) W, (35) 

we obtain for the corresponding RS functions 

,(«+!, 1-/3) 



w n (x; a,/3) = wq (x; a,/3) — sinh2x (n + a — + 1) — ^ — — ^ — — , (36) 



where 

w (x;a,j3) = -(a + l/2)cothx + (/3- l/2)tanhx. (37) 

As before, when n exceeds the value 2a, the corresponding Schrodinger equation enters in a disconjugacy regime. 
Using 

{p(a,/3) /., _ /n+a\ _ (<*+!)„ 
pKffl |V> _ r(2n+q+,g+l) n i Q ( „n-l\ ^ 8 i 

f ™ W ~ 2"n!r(n+a+,3 + l) X + U \ X )> 

from Eq (|33"|) and Eq (|34)) . we deduce for n > /3 — a — 1 

^(x;a,/3) ~ ^1/2(2+^ _> + 



^„(x;a,/3) - e (2n +Q -^+l )a; (2n ± «-^(n ± a-£ ± l) ^ +qo _ W 



i/> n having the same sign at both ends of the definition interval, ip n has no zero on R. The DBT built on the 
corresponding RS functions generate isospectral extensions of the HDPT potential. They are given by 

V^ [ x . a , 0) = V(x; a, 0) + w' n (x; a, 0) , (40) 

where w' n is given in Eq (|36| . Since l/ip n (x;a,/3) is divergent at the origin, it is not a physical eigenstate of V^ 
which is then strictly isospectral to V. It admits for bound state eigenfunctions 

tpf' (x;a,0) = (w n (x;a,0) - w k (x; a, /?)) ip k (x; a, 0) = ip o (x;a + l,0 - 1) % k _ 0) , (41) 

Pn ' (z) 

where k < a and where 

P^f (z) = (k + a-P + l) P,C-.-« (,) p£?'-W (z)-(n + a-P + l) P^_f^ +l) (z) PJ»-» (z) (42) 

is a polynomial of degree n + k — 1. 



B. Case of dispersion relation which is a second degree Laurent polynomial 

1. Eckart potential 
The Eckardt potential (a; G ]0,+oo[) can be written as jlL l2l. l28| 

V(x;a,b) = a(a + l)y 2 - 2by + V (a,b), a 2 < b, a,b > 0, (43) 

where y — cothx and 

V{a,b) = \+a. (44) 

The dispersion relation corresponds to a second degree Laurent polynomial of the form (a„ = a + n) 

9 b 2 1 b 2 
E n = a 2 + — - at - — (45) 

a < 

^-^(n + 2a) (n+ (a- - ) ) (n+ (a+- 



(n + a) z 
The bound states are obtained on the increasing part of E n , that is for < n < \[b — a and they are given by 

i/j n (x;a,b) = ip (x;a n ,p n )Pt^(y), (46) 

where 

^o(x;a,b) = {y-iy (y+1)* = e- faa; / a sinh a a;, (47) 

with 



&n ^*n 1 n 



Pn &n a 

The corresponding RS functions are 



(48) 



Pt- M (y) 

where 



„(x;a,6) =w (a;;an,Pn) + (2/ - 1J 5 „(a„.^„) ; ; — ' ^ 9 > 



a 

w o (x;a,0) = -(y + l) + ^(y-l). (50) 

Beyond the value n = a — Vo, the ip n do not correspond anymore to physical eigenstates and for n > - — a, E n 
becoming negative, the corresponding Schrodinger equation enters in a disconjugacy regime. 

Using Eq (j3"8)) , from Eq (|4*5)l we deduce the following asymptotic behaviour for i/) n (a„ + n = —2a — 2n) 

ip n (x;a,b) ~ (-2a)...(-n-2a+l)a;-( 2 " +a " +/3 ")/ 2 

II l\ I i\^ (a„ + l)...(a„+n) (51) 

ip n {x;a,b) ~ (y-1) 2 ^ — ^^ — ^ 
When 7?, > - — a, a n + n < and a n + /3 n = —2a — 2n. Consequently 



ip n (x; a, b) — > ±00 

X-S-0+ 

ip n {x;a,b) -> ± : 



± (52) 



with ± = (-1)". 

Since ip n has the same sign at both ends of the definition interval, the disconjugacy of the Schrodinger equation 
that it satisfies implies that ip n has no zero on R. We can then use the corresponding RS functions to build isospectral 
extensions of the Eckart potential given by 

7 W(g; a, b) = V(y; a,b)-2 (y 2 - l) *»» <* °» b) , (53) 

where w n is given in Eq (|49[) . l/ip n (x; a, b) is divergent at +00 and is not a physical eigenstate of V^ n > which is then 
strictly isospectral to V. The bound state eigenfunctions of V^- n ' are given by 

P (»f) ( z ) 
V'i™ ) {x; a, b) = (w n (x; a,b) - w k (x; a,b)) tp k (x; a,b) = Vo (x;a fe ,/3 fc ) — "; fc , (54) 

2PA " (2) 

where fc < v& — a and 

^ (*) ^ K (y + 1) + & (y 1)) P<£"» M (Z) Pt kA) (Z) (55) 

+ (y 2 - 1) (n + «„ + & + 1) P^ A) (2) P^»+ 1 ^+ 1 ) (z) 

- [fe # n] . 

2. Hyperbolic Rosen-Morse (HRM) potential 

Finally consider the HRM potential (x e R) which is given by [l], y, (2a] 

V(x;a,6) = a(a + l)y 2 + 2oy + V (a,o), a 2 > b, a,b > 0, (56) 

where t/ = tanh x and 

o 2 
V(o,6) = - 5 -o. (57) 

As for the Eckart potential, the dispersion relation corresponds to a second degree Laurent polynomial 

b 2 b 2 

E n = a 2 + — -a 2 n - — (58) 

a a n 

n , n J ( b\\ ( ( b 

— 2 \ n ~ 2a) In— a I I In— I a H — 



(n — a) 



(a n = a — n), but _E„ is now singular at the positive value n = a. 

As in the preceding cases, the bound states are obtained on the increasing part of E n , that is, for < n < a — yb. 
They are given by 

tl> n (x;a,b) = ^ fe^WP^' (y) , (59) 

where 



a e -bx/a 

TPo(x;a,b) = (l-y)^(l+y)* =—^-, (60) 



with 



O^n tin i ' 

6" 



Pn ^n a 

The corresponding RS functions are 



(61) 



( K\ ( a \ h *\ ^ + ^+(3n + lP i n a - 1 +1A + 1) (y) (R0 s 

{x;a,b) = w a (x;a n ,p n ) - (1 - y ) , (62) 

2 tf»A' (y) 



where 



Wo ( X ;a,f3) = ~(l + y)-^(l-y). (63) 

Beyond the value n = a — yb, the r/; n do not correspond anymore to physical eigenstates and -E„ is negative when 
n belongs to the intervals [a — b, a + b] and [2a, +oo]. 

using gum 

f p( Q <^) (i) = ( a+1 )^ 

1 Pt^i-i) = (-if^s (64) 

from Eq(|62j) we deduce the following asymptotic behaviour for ip n 



X— ¥ — OO 



4>n(x;a,b) ~ (1 _ y) ^ (a„ + l) (a„+n) _ ( 65 ) 



J' — >■ + OO 



Consider first the case (i) where a — - < n < a. Then a n > 0, a n > and j3 n + n < 0. Consequently 

( ?/; n (#; a, 6) — > -hoo 



a;— )•— oo 



) *P»(x;a,b) -► 0+. (66) 

^ x— >+oo 

Consider now the case (ii) where a + - > n > a (a n < 0). In this case we have a n + n < and /3 n > 0. It results 

( \j) n (x; a, 6) — > ±oo 



x— J- — OO 

ip n (x;a,b) -> 0*, 

a;— J-+oo 



± (67) 



where ± = (—1)". 

Finally consider the case (hi) n > 2a (<z„ < 0). We have a n , Pn < and a n + n, j3 n + n > 0. ip n is then divergent 
both at +oo and — oo. Depending on the value of n (compared to a + b/k, k 6 N*), we can then have the same sign 
or not for the limits of ip n . 

When tp n has the same sign at both ends of the definition interval, the disconjugacy of the Schrodinger equation 
that it satisfies implies that ip n has no zero on K. We can then use the corresponding RS functions to build isospectral 
extensions of the HRM potential given by 

V (B) (x; a,b) = V(y;a,b) + 2(l- y 2 ) dWn ( . y; °' &) , (68) 

ay 



10 

where w n is given in Eq (|62l) . 

In the cases (i) and (ii), l/ip n (x;a,b) is divergent at +00 or — 00 and is not a physical eigenstate of V' 71 ' which is 
then strictly isospectral to V. 

At the contrary in the case (iii), when l/if) n (x;a,b) is regular, it also satisfies the required Dirichlct boundary 
conditions and constitutes the fundamental bound state of V^ . The DBT is then a backward SUSY partnership. 
We will not consider this case in the following. 

In the cases (i) and (ii), the bound state eigenfunctions of V"( n ) are given by 

ip^\x;a,b) = (w n (x;a,b) - w k (x;a,b))ip k (x;a,b) = V>o (x]a kl f3 k ) "; ' . , (69) 

2P y n n ' Pn > (z) 



where k < a — \fb and 



P { :f (*) = (1 - y 2 ) (k + a k + h + 1) PMJ (z) pjp+^+Q {z) (70) 

- K (i + v) - fa (1 - v)) PL an ' M (*) P ( k akA) (*) 

- [k O n] . 



IV. SHAPE INVARIANCE 

For the isotonic oscillator we have proven explicitely in [25l . |27| that the shape invariance property of the initial 
potential is transmitted to all its strictly isospectral (LI and L2 series) successive extensions. This is not the case for 
the extensions obtained from the potentials considered in [26| via the use of regularizing symmetries. 

To look for such a property for the extensions Morse and HDPT potentials obtained above, we have to consider the 
superpartner of the potential V^ (x; a) which is defined as 

V' n » (s; a) = V< n > (x; a) + 2 (w^ (x; a)) ' , (71) 

he RS function as 
Since (see Eq©) 



Wq (x; w, a) being the RS function associated to the ground level of V^ 



(n), x / n E n (a) 

w o \x;a) = -w n (x;a) r r, 

w n (x;a) - wo(x;a) 



we have, using Eq([T6 



V^(x;a) = V(x;a)-2(—— 
\w n [x; 



E n {a) 



(72) 



(73) 



;o) -w (x;a) 
We suppose that V is a TSIP satisfying 

V(x; a) = V(x; a) + 2w' {x; a) = V(x; a - 1) + Ei(a). (74) 

Inserting Eq([?l|) into Ea([73]). it results 

V in) {x;a) = V(x;a-l) + E 1 (a)-2W'(x;a), (75) 

where 

W(x; a) = w (x; a) H r^ -. (76) 

w n [x\a) -w (x;a) 
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A. Morse potential 



Consider the case of the Morse potential. Using Eq (|22)) and Eq ([23|) . Eq (f76| gives (z = 2by) 



Ll a .\(z) 



W(x;a) = -a n + by + z ^ // (77) 



and since a n = (a — l) n _ 1 , we deduce 



W(x; a) = - | w a (x; (a - 1)„_ 1 ) - z " 2(a „ 1)ni j = : -w n ^{x; (a - l)^). (78) 

Ai-1 " 1 ( z ) 



Consequently 



V< n >(a;;a) = y(a;;a-l)+2<_ 1 (.T;(a-l) n _ 1 )+i; 1 ( a ) 

= V^-^ix-a^+E^a). (79) 

This is not strictly speaking a shape invariance in the sense of Gendenshtein [3|. As noted by Quesne [29(, we 
rather obtain a kind of enlarged shape invariance property where the SUSY QM partner of the n th extended potential 
V^ n ' has not the functional form of V^- n ' (with translated parameters and an additional constant) but the one of the 
preceding extension U^™ -1 '. 

B. HDPT potential 

For the HDPT potential, combining Eq (|36|) and Eq (|37|) in Eq (l76|) . we find that W takes the form (z = cosh 2x) 

AnP( a, ~^ (z) 

W(x;a,P)=w (x;a,P)+ ■ " - , K _> . (80) 



(a + l,l-/3) 



smh2xP^_T''- p ' (z) 



But P„ ' (z) satisfies the differential equation 



(1 - z 2 ) y"{z) + ( a + j3-z(a-!3 + 2)) y'{z) + n (n + a-0 + l) y{z) = 0, (81) 

which, combined to Eq ([551) gives 

4nPi a -V (z) = (z 2 - 1) (n + a - /? + 2) p(«+ 2 <^+ 2 ) (z) + 2 (a + P - z (a ~ P + 2)) p^f'-P+V ( z ) . (82) 

Substituting this result in Eq ([8"0)) and using Eq (l571) . we obtain 

W{x-a,P) = ^— + sinh 2a: (n + a-ff + 2) {a+1A _ p) (83) 

= -w n _i (»;ai,/3i), 
that is, 

U<">(z;a,/3) = V(a;;oi, i 9i) + 2«;;_ 1 (a;;ai >i 8i) + £;i(o) 

= ^"-^(ajjai^O+E^o). (84) 

Again, yw satisfies the enlarged shape invariance property defined above. 
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V. CONCLUSION 

For PTSIP with a finite number of bound states, we have enlarged the generalized SUSY QM partnership presented 
in [2414261 128[ by showing another way to obtain disconjugated unphysical eigenfunctions which may serve to generate 
regular rational extensions of these potentials. We have studied more explicitely those of these extended potentials 
which are strictly isospectral to the initial one. These results encompass in particular those obtained very recently 
by Quesne [29[ for the Morse potential. The enlarged shape invariance property revealed by Quesne for the Morse 
potential extensions is proven to be shared by the HDPT potential. 
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